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ABSTRACT 

The recent observations of cross temperature-polarization power spectra of the Cosmic 
Microwave Background (CMB) made by the WMAP satellite are in better agreement with a 
high value of the Thomson scattering optical depth r w 0.17. This value is close to r = 0.3, 
which is taken as the upper limit in the parameter extraction analysis made by the WMAP 
team. However, models with r ~ 0.3 provide a good fit to current CMB data and are not 
significantly excluded when combined with Large Scale Structure data. By making use of 
a self-consistent reionization model, we verify the astrophysical feasibility of models with 
t ~ 0.3. It turns out that current data on various observations related to the thermal and 
ionization history of the intergalactic medium are not able to rule out r ~ 0.3. The possibility 
of a very extended reionization epoch can significantly undermine the WMAP constraints on 
crucial cosmological parameters such as the Hubble constant, the spectral index of primordial 
fluctuations and the amplitude of dark matter clustering. 



1 INTRODUCTION 

The recent results on Cosmic Microwave Background (CMB) 
anisotropy from the WMAP satellite (Bennett et al. 2003; Spergel 
et al. 2003) represent an extraordinary success for the standard cos- 
mological model of structure formation based on Cold Dark Matter 
(CDM) and adiabatic primordial perturbations. Furthermore, data 
releases from the Sloan Digital Sky Survey (SDSS; Tegmark et al. 
2004) are living up to expectations and combined analysis of all 
these datasets are placing strong constraints on most cosmological 
parameters. 

However, even if theory and observations are in spectacular 
agreement, discrepancies seem to be present and have already stim- 
ulated the interest of several authors. One of most unexpected re- 
sults from the first year WMAP data is indeed the detection of 
an excess in the large angular scale cross temperature-polarization 
power spectrum (Kogut et al. 2003). In the standard scenario, this 
excess may be interpreted as due to a late reionization process 
which produces a second last scattering surface at lower redshifts 
(~ 10 — 20). The CMB radiation can therefore get polarized 
through anisotropic Thomson scattering on angular scales compa- 
rable to the horizon at reionization (0 > 5°). 

Analysis made by the WMAP team of the cross temperature- 
polarization data constrains the best-fit value of the Thomson opti- 
cal depth of the universe to r ~ 0.17 (the exact best-fit value and 
the errors depend on the analysis technique employed). Following 
this, considerable effort is spent in constructing models, both nu- 
merical (Ciardi, Ferrara, & White 2003; Ricotti & Ostriker 2004; 
Sokasian et al. 2003; Sokasian et al. 2004; Gnedin 2004) and semi- 
analytical (Wyithe & Loeb 2003; Hui & Haiman 2003; Choud- 
hury & Ferrara 2005), analyzing the feasibility of such a high r. 
In this context, it was shown that an early population of metal-free 
stars, with reasonably low values of star-forming efficiency and es- 



cape fraction (Ciardi, Ferrara, & White 2003; Choudhury & Ferrara 
2005), is sufficient to produce such a high r. 

In this paper we want to emphasize another aspect of this 
result and investigate the compatibility of cosmological and as- 
trophysical data with a very extended reionization process i.e. 
r ~ 0.3. As we illustrate in the next section, the parameter extrac- 
tion made by the WMAP team has not considered this case, since 
a top-hat prior of r < 0.3 has been included in their analyses (see 
Spergel et al. 2003 and Peiris et al. 2003). Given the fact that the 
reionization process is still to be fully understood, one could won- 
der about the robustness of this assumption. In the next section we 
will indeed show that r = 0.3 is compatible with the data and, 
as already pointed out in (Tegmark et al. 2004; Fosalba & Szapudi 
1998), letting r to vary freely up to r = 0.5 more than doubles 
the error bars on cosmological parameters. Studying the astrophys- 
ical feasibility of a very extended reionization is therefore not only 
important for the theoretical understanding of the process but also 
for determining the cosmological parameters with highest possible 
accuracy and reliability. 

The paper is organized as follows: In the next section, we ana- 
lyze current CMB and Large Scale Structure data to show that high 
values of the optical depth are not excluded and discuss the impli- 
cations of having r ~ 0.3 on cosmological parameters. In section 
3, we verify whether models with r ~ 0.3 violate any observa- 
tional constraints related to the ionization and thermal history of 
the universe. The final section presents our conclusions. 



2 COSMOLOGICAL CONSTRAINTS ON THE OPTICAL 
DEPTH 

In this Section we investigate the possibility of having a high op- 
tical depth r ~ 0.3 by analyzing current CMB and Large Scale 
Structure data. The method we adopt is based on the Markov Chain 
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Figure 1. Likelihood distribution functions on the cosmological parameters described in the text. Only the WMAP data is used without prior on the optical 
depth. As we can see, values of r ~ 0.3 are well consistent with the data 
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Figure 2. 1-<t and 2-cr likelihood contour plot in the Hq vs r plane from 
the WMAP data alone. There is a clear degeneracy between the two param- 
eters and the Hubble parameter is unconstrained from the CMB data alone. 
Including the prior r < 0.3 artificially constrains the Hubble parameter. 



Monte Carlo package cosmomc 1 (Lewis & Bridle 2002). We use 
individual chains and we adopt the Raftery and Lewis convergence 
diagnostics 2 We sample the following 6 dimensional set of cos- 
mological parameters, adopting flat priors on them: the physical 
baryon and CDM densities, u>i, = Qth 2 and tJdm = SldmA 2 , the 
ratio of the sound horizon to the angular diameter distance at de- 
coupling, 8 S , the scalar spectral index and the overall normaliza- 
tion of the spectrum, n s and A s , and, finally, the optical depth to 
reionization, r. Furthermore, we consider purely adiabatic initial 
conditions, we impose flatness and we do not include gravitational 
waves. Relaxing flatness does not affect r in a very significant way; 
however, it does alter the values of the other basic parameters and 
thus, may further change the interpretation and astrophysical im- 
plications of the overall results. We restrict our analysis to the case 



1 Available from http : //cosmologist . info. 
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of 3 massless neutrino families. We include the first-year tempera- 
ture and polarization data (Bennett et al. 2003) with the routine for 
computing the likelihood supplied by the WMAP team (Verde et 
al. 2003) as well as the CBI (Readhead et al. 2003) and ACBAR 
measurements of the CMB. 

In addition to the CMB data, we also consider the constraints 
on the real-space power spectrum of galaxies from the SDSS 
(Tegmark et al. 2004). We restrict the analysis to a range of scales 
over which the fluctuations are assumed to be in the linear regime 
(k < 0.2/i _1 Mpc). When combining the matter power spectrum 
with CMB data, we marginalize over a bias b, constant with scale 
and redshift, considered as an additional nuisance parameter. 

The results of our analysis are reported in Table 1 and dis- 
played in Figures 1^1. In Figure 1 we show the likelihood distribu- 
tion function obtained by analyzing the WMAP data only. As we 
can see the data does not provide strong constraints on r and values 
as large as r ~ 0.5 are compatible with the data. As a further ex- 
ample, we plot in Figure 2 the 2 — D likelihood distribution on the 
Hq-t plane after marginalization over the remaining nuisance pa- 
rameters. As we can see there is a clear degeneracy between the two 
parameters and the value of the Hubble constant is not constrained 
by the WMAP data alone. The claimed constraint in Spergel et al. 
(2003) of h = 0.72 ± 0.05 at 68% c.l. is mostly due to the top prior 
on r. For comparison, we report in Table 1 the results from a simi- 
lar analysis but with now the inclusion of an external top-hat prior 
of t < 0.3 as in the analysis made by the WMAP team (Spergel 
et al. 2003). As we can see, the inclusion of the optical depth prior 
more than halves the error bars and greatly improves the constraints 
on the Hubble parameter, the spectral index and the baryon density 
to values in agreement with those reported by the WMAP team. 

It is useful to compare the WMAP constraints derived in the 
case of t ~ 0.30 with independent cosmological probes. If we 
assume a high value of r = 0.30 ± 0.01, then the likelihood 
best fit values are u>t = 0.026, h — 0.78 and n s — 1.07. A 
value of h ~ 0.8 is still in reasonable agreement with the HST 
constraint of h = 0.72 ± 0.07 (Freedman et al. 2001). Values of 
cut, > 0.025 are difficult to reconcile with abundances of primor- 
dial light elements and standard Big Bang Nucleosynthesis which 
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provide oj b = 0.020±0.001 (Buries & Tytler 1998). However, neu- 
trino physics is still rather unknown and deviations from the stan- 
dard model, like inclusion of neutrino chemical potential or extra 
relativistic degrees of freedom could be present and larger values 
of u>b may be possible (see e.g. Hansen et al. 2002). We can there- 
fore conclude that the high r case is in agreement with the WMAP 
data and that the best fit cosmological model with r = 0.3 is in 
reasonable agreement with several other observables. 

Including the SDSS data, as we do in Figure 3, breaks several 
degeneracies, mainly between the spectral index n s and r, and im- 
proves the constraints. However, as we can see from Figure 4, mod- 
els with r = 0.3 are still inside the 2-a level and are not strongly 
disfavoured. In the WMAP+SDSS analysis, the best fit model un- 
der the assumption of r = 0.30 ± 0.01 has u> b = 0.027, h = 0.75, 
n 3 = 1.08 and as = 1.17. These results are close to those from the 
previous WMAP analysis further confirming the reasonable agree- 
ment between the SDSS data and r = 0.3. 

Including ACBAR and CBI experiments to WMAP+SDSS 
improves the constraints on r and we get an upper limit of r ^ 
0.25. While this disfavours the r ~ 0.3 hypothesis, one should 
also note that other effects like massive neutrinos (see Fogli et al. 
2004), running spectral index and a scale-dependent bias, just to 
name a few, could all be present and improve the agreement be- 
tween models with high r and the data. For example, including a 
running spectral index in the primordial spectrum, which may be 
motivated by inflation and/or by a scale dependent bias in galaxy 
clustering, enlarges the constraints again to r < 0.32. In the case 
of running spectral index and r = 0.30 the best fit parameters are 
uj b = 0.024, h = 0.67, n s = 0.93, dn/dlnfc = -0.08 and 
a 8 = 1.11. 

An accurate determination of as may be the best way to deter- 
mine r. Unfortunately there is a large scattering in the constraints 
on as from current observations. A value of erg = 1.1 — 1.2 is 
compatible with the determinations inferred from the local clus- 
ter X-ray temperature function (see e.g. Pierpaoli et al. 2001; Bor- 
gani & Guzzo 2001; Eke et al. 1996) and cosmic shear data (Van 
Waerbeke et al. 2001; Refregier et al. 2002; Bacon et al. 2003), 
but incompatible with other analysis (see e.g. Seljak 2002; Viana 
et al. 2002). Recent measurements of galaxy clustering at redshift 
z = 0.6 from the COMBO-17 survey gives cr 8 = 1.02 ± 0.17 at 
zero redshift (Phleps et al. 2005). It is interesting to point out that, 
if the small-scale CMB anisotropy excess measured by the CBI 
experiment (see Readhead et al. 2003) is indeed due to integrated 
Sunyaev-Zeldovich effect, this would be compatible with a value of 
as = 1. 04±0. 12 (Komatsu & Seljak 2002) therefore in agreement 
with a very extended reionization scenario. 



3 ASTROPHYSICAL CONSTRAINTS 

A potential challenge to models with high values of r can come 
from the astrophysical constraints. To address this issue, we study 
the implications of a high value of r ~ 0.3 on reionization his- 
tory of the universe. For this purpose, we use a semi-analytical 
model (Choudhury & Ferrara 2005) developed for studying cos- 
mic reionization and thermal history of the intergalactic medium 
(IGM). The model implements most of the relevant physics govern- 
ing these processes, such as the inhomogeneous IGM density distri- 
bution, three different classes of ionizing photon sources (massive 
PopIII stars, PopII stars and QSOs), and radiative feedback inhibit- 
ing star formation in low-mass galaxies. The main advantage of the 
model is that it can be constrained quite well by comparing its pre- 



dictions with various observational data, namely, the redshift evo- 
lution of Lyman-limit absorption systems (Storrie-Lombardi et al. 
1994), Gunn-Peterson (Songaila 2004) and electron scattering op- 
tical depths (Kogut et al. 2003) and cosmic star formation history 
(Nagamine et al. 2004); in addition the model also includes con- 
straints obtained by Salvaterra & Ferrara (2003) on Near Infra Red 
Background (NIRB) data (Matsumoto et al. 2000; Gorjian, Wright, 
& Chary 2000; Cambresy et al. 2001; Wright 2001). In this pa- 
per, we extend the above model of Choudhury & Ferrara (2005) 
to include the possibility of a high r along with high values of 
n s — 1.08, as = 1.18 and ojt, = 0.026, as discussed in the previ- 
ous section and reported in Table 1 . 

The most crucial difficulties in admitting such a high value of 
t within standard reionization scenarios can be divided into two 
parts. In the first part, we discuss whether it is possible to have 
sources which can produce a high number of ionizing photons so as 
to produce a large r. A common mechanism for obtaining a large 
t is through metal-free PopIII stars, which has been used widely 
(see, for example, Cen 2003; Cen 2003; Venkatesan, Tumlinson, 
& Shull 2003), and is in agreement with various set of observa- 
tions (Choudhury & Ferrara 2005) including the NIRB data (Sal- 
vaterra & Ferrara 2003; Magliocchetti, Salvaterra, & Ferrara 2003). 
In such scenarios, PopIII stars with a photon production efficiency 
em = esF,ni/esc,in ~ 0.006 (where esF,m is the star-forming 
efficiency and / csc ,iii is the escape fraction of ionizing photons) 
are able to produce a r = 0.17 (the best-fit value obtained by fit- 
ting the temperature-polarization cross power spectrum to ACDM 
models in which all parameters except r assume their best fit values 
based on the temperature power spectrum; Kogut et al. 2003). One 
might then expect that a considerable higher value of em must be 
required to produce a r of 0.3. However, one should also keep in 
mind that the model with r = 0.3 has higher values of erg = 1.18 
and u)b = 0.026 compared to the standard one, which means that 
(i) the number of collapsed haloes is higher and (ii) the baryonic 
mass within the collapsed haloes is larger. In fact, the value of em 
required to match the high r is w 0.0035, which is similar to the 
value in standard scenarios. So we conclude that production of a 
large number of ionizing photons at high redshift does not pose any 
serious difficulty for models with high r. 

Various other mechanisms for early reionization have been 
suggested in order to produce a high r, most notably being an early 
population of accreting black holes (Madau et al. 2004; Dijkstra, 
Haiman, & Loeb 2003; Ricotti & Ostriker 2004). The accretion of 
gas onto these black holes or miniquasars would produce a back- 
ground of hard photons, which have a large mean free path and can 
ionize large regions of the IGM. However, the constraints from the 
present-day soft X-ray background imply that these sources can, at 
most, produce ~ 3 photons per hydrogen atom (Dijkstra, Haiman, 
& Loeb 2003; Salvaterra, Haardt & Ferrara 2005), which is much 
less than the requirement for fully ionizing the IGM - hence the 
reionization from these sourced would be partial. While the nature 
of reionization (topology of the ionized regions, thermal state of 
the IGM etc) would be quite different if it is achieved by the black 
holes rather than PopIII stars, whether the photons from the mini- 
quasars can produce a high value of r ~ 0.3 without violating any 
other observational constraints needs to be verified. 

The second difficulty regarding high values of r is whether 
such a high number of photons at high redshift violate any obser- 
vational constraints (see Figure 5). The first set of such constraints 
come from the observations related to the IGM. Since there are vir- 
tually no observational data on IGM at z > 6, the most severe 
constraints on the model comes from the Lyman-a Gunn-Peterson 
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Figure 3. Likelihood distribution functions on the cosmological parameters described in the text. The WMAP+SDSS data is used without external prior on 
the optical depth. Values of r ~ 0.3 are now less in agreement with the data but still inside the 2cr level. 
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Figure 4. 1-cr and 2-cr likelihood contour plots in the r vs erg (left panel) and r vs n s (right panel) planes from a WMAP+SDSS analysis. Even after the 
inclusion of large scale structure data, there is a clear degeneracy between the parameters. Values of r ~ 0.3 are in between the 2 a contours if n s ~ 1.1 and 
cr 8 ~ 1.17. 



Table 1. Mean values and marginalized 95% c.l. limits for several cosmological parameters from WMAP and SDSS (see text for details). 



Parameter 


WMAP only 


WMAP + r < 0.30 


WMAP + r = 0.30 ± 0.01 
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WMAP + SDSS + t = 0.30 ± 0.01 
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optical depth tgp at z « 6. In the standard scenario, tgp obser- 
vations are very well matched by the flux from normal PopII stars 
with an efficiency en « 0.0005 along with a population of QSOs 
as required by data on optical luminosity function. It is thus nec- 
essary that the PopIII stars do not have a significant flux at 2 ~ 6. 
For the models studied in this paper, this provides a mild constraint 
on the transition redshift of PopIII stars Ztrans > 10.5. However, 
there is another physical process which affects tgp and thus pro- 
vides further constraints. This has to do with the mean free path 
Ah of hydrogen-ionizing photons. At high redshifts, the high flux 
from PopIII stars create huge ionized regions, thus producing high 
values of Ah ~ 60-80 Mpc. Once the PopIII stars disappear at 
Ztrans, these ionized regions start becoming neutral over a recom- 
bination time-scale resulting in a gradual decrease of Ah- However, 
the constraints on tgp (taking into account the uncertainty in the 
value of the slope of the IGM temperature-density relation 7) im- 
ply that Ah < few Mpcs at 2 ~ 6. To obtain such low values of 
Ah, one has to make sure that the PopIII stars start disappearing 
around 2 tra ns ~ 11- Once this condition is satisfied, the model is 
consistent with various available observations, namely, the redshift 
evolution of Lyman-limit systems, the temperature corresponding 
to the mean IGM density, the cosmic SFR and tgp, as shown in 
Figure 5. 

The second set of observational consequence is related to the 
NIRB which is believed to be due to the PopIII stars. The obser- 
vations of NIRB constrain the combination esF,m(l — jcsc.in) 
and ztrans (Salvaterra & Ferrara 2003; Magliocchetti, Salvaterra, 
& Ferrara 2003). First, note that the value of r is sensitive only to 
the combination esF,in/esc,ni; hence, it is not difficult to match the 
NIRB constraints on esF,m(l — /csc.iii) by suitably choosing the 
values of esF.in and / csc ,iii. Second, the decline of the NIRB at 
wavelengths below the J band constrains 2 tr ans ~ 10, whereas the 
tgp constraints at z > 6 require that ztrans > 11 (as discussed in 
the previous paragraph). However, the constraint from NIRB is not 
too strong mainly due to lack of reliable data at wavelengths short- 
wards of the J band. It is thus difficult to rule out models where the 
PopIII star formation rate starts decreasing gradually from 2 ~ 11 
with substantial contribution remaining till 2 ~ 9. 

In conclusion, current observations, related to the IGM and 
NIRB, are not able to rule out models with r ~ 0.3, though they 
are only marginally consistent in some cases. 



4 CONCLUSIONS 

In this paper we have investigated the feasibility of a very extended 
reionization. We can summarize our results in 

• The first year WMAP data is in perfect agreement with an 
Thomson scattering optical depth parameter as high as r = 0.5. 
Not considering this possibility in parameter extraction analysis 
may bias the results on several quantities like the Hubble constant, 
the baryon density and erg. 

• Including galaxy clustering data from SDSS disfavours mod- 
els with r = 0.3 but not at high significance. Similarly, inclusion 
of the CMB data from CBI and ACBAR improves the constraints 
tor ^ 0.25. On the other hand, if additional parameters like a run- 
ning spectral index or scale dependent bias are considered, high r 
models are found to be in better agreement with observations. 

• We have therefore studied the astrophysical feasibility of such 
a high value of r, particularly with respect to the observations re- 
lated to the reionization history of the universe. Using a simple 
semi-analytical model (Choudhury & Ferrara 2005), we find that 
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Figure 5. Comparison of a model having r = 0.3 with various as- 
trophysical observations. Adopted parameters are em = 3.5 X 10~ 3 , 
en = 3 X 10~ 4 , ztrans = 11. The panels show as a function of redshift: 
(a) filling factor of ionized hydrogen and doubly-ionized helium regions, (b) 
specific number of Lyman-limit systems, (c) ionizing photons mean free 
path for hydrogen and helium, (d) Lyman-ct Gunn-Peterson optical depth 
for three values of slope 7 of the IGM temperature-density relation. 

an early population of massive metal free PopIII stars with photon 
production efficiency em = esF,ni/csc,iii ~ 0.0035 is sufficient 
to produce such a high r. However, in order that the photon flux 
does not violate the Lyman-a Gunn-Peterson optical depth con- 
straints at 2 > 6, the PopIII star formation rate should start de- 
creasing around Ztrans ~ 11- This value of 2 tra ns is marginally 
consistent with the observations of NIRB. 

We can therefore conclude that the possibility of a very ex- 
tended reionization epoch is still not completely ruled out by cur- 
rent data and should be taken into account, especially when deriv- 
ing reliable 2a constraints from CMB anisotropy measurements. 
In future, more accurate CMB polarization measurements, such as 
those expected from the Planck satellite experiment, together with 
high accuracy measurements of erg from weak lensing and galaxy 
surveys, will be able to verify or rule out this interesting possibility. 
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